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Abstract. 

We study ultracold atom-ion collisions in the presence of an external magnetic 
field. At low collision energy the field can drastically modify the translational motion 
of the ion, which follows quantized cyclotron orbits. We present a rigorous theoretical 
approach for the calculation of quantum scattering amplitudes in these conditions. 
Collisions in different magnetic field regimes, identified by the size of the cyclotron 
radius with respect to the range of the interaction potential, are investigated. Our 
results are important in cases where use of a magnetic field to control the atom-ion 
collision dynamics is envisioned. 
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1. Introduction 

The central role played by the coUisional interaction in determining the macroscopic 
properties of cold atomic gases has motivated extensive theoretical and experimental 
studies of neutral atom collisions at very low collision energies over the last decade. For 
instance, magnetic Feshbach resonances have been discovered and theoretically analyzed 
in several atomic species (see, e.g., [D l2])- Such resonances allow the interatomic 
interaction to be finely tuned, and have been used to explore a variety of phenomena 
such as the BEC-BCS crossover [3J and strongly correlated quantum phases 

While atom-ion collisions at standard temperatures are a well established subject 
[5j, their behaviour at very low collision energy has only recently begun to raise a rapidly 
increasing interest. Atom- ion systems are in particular very promising candidates for 
the realization of coUisional quantum gates [6] and for the study of ionic impurities 
in neutral Bose- Einstein condensates [7]. Experimental progress makes now ultracold 
atom- ion systems addressable in a controlled way [8j. First measurements of elastic 
cross sections and sympathetic cooling rates between a Bose condensed gas and a single 
ion held in a Paul trap have been reported in Refs. [8], [9] . Charge exchange has also 
been theoretically analyzed both for its interest in fundamental molecular physics and 
in order to quantify harmful sources of atom-loss or decoherence in the experiments 

mm- 

The presence of magnetically tuned Feshbach resonances induced by the hyperfine 
interaction has been demonstrated in free space in [10]. Such resonances could in 
principle be engineered to increase the rapidity of coUisional quantum gates or to 
manipulate the many-body properties of a neutral gas in the presence of a charged 
impurity [12]. However, non trivial effects have to be taken into account to describe 
quantitatively atom-ion collisions in most non idealized situations, such as trap 
confinement [13]. An additional feature neglected in previous work [10] is the influence 
of the magnetic field on the translational motion of the ion, which will classically follow 
helicoidal rather than rectilinear trajectories. Quantum mechanically, motion in the 
plane transverse to the field has to be quantized, resulting in the discretization of the 
continuum in Landau levels [H]. 

The dramatic influence of even a weak magnetic field on electron photodetachment 
was observed a long time ago [15j and modeled using a frame transformation 
approximation, which relies on the much smaller length scale of the interatomic potential 
compared to the cyclotron radius fT6]. However, the frame transformation approach 
cannot be applied to ultracold atom-ion collisions since the typical length scales become 
comparable and cannot therefore be separated already at commonly used field intensities 
of a few hundred Gauss (see below). Numerically exact i?- matrix calculations have been 
carried out for hydrogen photoionization but were limited to strong magnetic fields |17j . 
The main goal of our paper is to develop a rigorous theory applicable to atom-ion 
collisions in a magnetic fields of arbitrary strength. 

The paper is organized as follows. The system and general scattering properties in 
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a magnetic field are introduced in Sec. 2. Sec. 3 presents the numerical methods. In 
Sec. 4 our framework is applied to the scattering of Ca+ ions and Na atoms. A variety of 
threshold, quantum interference and resonance phenomena are identified and discussed. 
Concluding remarks and perspectives end this work. 



2. Theoretical model 

2.1. Background 

We consider a low energy collision between a neutral atom with a singly charged ion 
in the presence of a uniform static magnetic field B oriented in the z direction. We 
suppose that both particles have spherical symmetry, vanishing internal spin, and that 
the atom has an infinite mass. The latter assumption allows the number of spatial 
degrees of freedom to be reduced from six to three, thus avoiding complications related 
to the separation of the center of mass motion [TSj . 

The spherically symmetric interparticle molecular potential V{r) has a long-range 
behaviour determined by the second order charge induced-dipole interaction — C4r~^, 
where 2C4 is the static dipolar polarizability of the neutral atom. The detailed shape 
of the potential well is immaterial for ultracold collisions and it can be to excellent 
approximation described by energy-independent quantum defects or scattering lengths 
[T3] . For this work we take therefore a simple analytical form V{r) = D^. {x~^ — 2x~^) 
with X = r/ve. The potential parameters De and are adjusted to approximately 
agree with the pubhshed CaNa"*" potentials and to reproduce the C4 = 81.3 a.u. value 
[llj . The s-wave scattering length is currently unknown for any atom- ion system. We 
perform therefore some exploratory studies for different scattering lengths by adding 
a small correction to the bottom of the potential well to vary the value of in the 
numerical calculations [T9] . 

The Hamiltonian for a singly positively charged ion is 

(p- eX) 

H-^^^V ^ (1) 

where A is the vector potential associated to B, p the generalized momentum and e > 
the ion charge. The cylindrical symmetry about the direction of the magnetic field 
implies conservation of the axial projection of the angular momentum L = f x p, 
whose quantized values will be noted ^A. It is convenient to choose the origin of the 
coordinates on the neutral atom and to work in the symmetric Landau gauge, in which 
the vector potential has the form A = ^B x r. With this choice, after some simple 
algebra the Hamiltonian takes in cylindrical coordinates {p, ip, z} the following form 



pI pi Ll muj"^ 

— + — H \ 

2m 2m 2mp2 2 



H=^ + -^ + + -ojh, + V (2) 



of a particle in a diamagnetic potential of Larmor angular frequency ^ = ^ interacting 
with a fixed center of force V . The coordinate p is the distance from the z axis, Pp is 
its conjugate momentum and Pz is the linear momentum in the direction of B. 
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Conservation of the axial angular momentum reduces the orbital Zeeman term wL^ 
to a constant but its influence on the spectrum of H is essential. In fact, quantized 
energies for ion motion in the plane orthogonal to the Landau levels [I^, can be 
expressed in terms of A and of the radial quantum number n by 

EnK = nw{2n+\K\- K + l) , n = 0,l,2,---. (3) 

Eigenvalues have therefore infinite degeneracy, a quantum counterpart to the classical 
freedom to choose the cyclotron orbit center keeping the energy fixed. The level spacing 
2hw for n — 7- nil transitions corresponds to the energy of the classical cyclotron motion 
with angular frequency twice the Larmor frequency of precession, Uc = eB/m = 2u . 

The explicit expression of the eigenfunctions for planar motion with azimuthal angle 
If is well known (see, e.g., Ref. [H]) 



where Oho = \ltifmhj is the harmonic oscillator length and Nnk a normalization 
constant. The quantum description of helicoidal motion is finally completed by a uniform 
translation in the z direction 

<^>fA(r) = 0nA(p,^)|A;„r^/V'="^ (5) 

Note that to simplify the notation we have defined fc„ = /cnA dropping the conserved 
quantum number A. Conservation of the total energy E determines the magnitude of 
the longitudinal wave vector 

E=^ + i?„A. (6) 

where the first term on the rhs is the kinetic energy in channel n, ii^kin = 2ni ' 

Classically, as the atom and the ion approach, the spherically symmetric interaction 
distorts the cyclotron orbit modifying the energy distribution between longitudinal 
and transverse degrees of freedom. Quantum mechanically, a Landau collisional state 
can be either refiected or transmitted across the collision center with different 
probability amplitudes. In general, a collision starting with initial quantum number 
n will populate all states ^^'a energetically accessible at energy E with a transition 
probability amplitude These energetically open channels are characterized by 

kf^i > 0, see Eq. ([6]). The closed channels k"^, < are not energetically allowed at large 
atom-ion separation but can be populated during the collision giving rise to dramatic 
resonance effects. 

Inspection of the typical length scales in the problem shows the essential role played 
by Landau quantization for cold atom-ion collisions. In the ultracold regime, the range 
of the interaction can be defined in terms of the polarizability coefficient C4 through 
R* = ^j2fmU^~^ [To]. The R* parameter is proportional to the location i?*/7r of the 
last node in the zero-energy wavefunction for infinite scattering length [19] . As such, it 
gives an estimate of the distance beyond which the wave function oscillates freely. The 
characteristic energy E* = / {2mR*'^) sets the scale at which quantum effects become 
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Table 1. Characteristic energy and length scales as defined in the text for few sample 
atom-ion systems in a magnetic field B — 500 G. The oscillator length at this value of 
the field is aho — 3066ao. Calculations use the reduced mass of the atom-ion dimer. 



System 


i?*(ao) nujc{,li^) E*{fiK) 


Na+Ca+ 
Rb+Ba+ 
Yb+Yb+ 


2081 2.5 1.4 
4765 0.70 0.92 
5572 0.43 1.1 



important [20]. The quantum length scale associated to cyclotron motion is Oho which 
corresponds to the size of the lowest energy Landau state. The Landau levels spacing 
is finally given by huc- 

Let us now consider a magnetic field intensity B = 500 G. Tab. [U shows the 
characteristic quantities for some systems of current or possible experimental interest 
at this value of magnetic field. One can remark that the range of the interaction is 
comparable to the harmonic oscillator length. Moreover, the transition between the 
weak R* ^ and the strong field R* ^ a^o limits can be explored by varying B over 
an experimentally accessible range. The Landau splitting of the order of fiK makes the 
discretized nature of the continuum in the magnetic field essential if collision energies 
of this order of magnitude are considered. 

The scattering wavefunction presents the asymptotic |2;| — > oo behaviour 

^nAi^ -> ^nAi^ + Yl (sign(^)A:n' ^ K) ^f^A (P, 1^1) (7) 

fc„,>0 

in terms of an incoming <l>^^ Landau state and of scattered waves outgoing in the 
open channels. One should note that since the transverse motion is confined by the 
magnetic field, scattering is effectively one dimensional. The measurable refiection 
and transmission coefficients for a progressive kn > incoming wave can be promptly 
identified from Eq. ([7]): 

Rn'n = l/n'ra (~^n' ^n) | ; ^n'n = \^n'n + fn'n (^n' ^n) | • (8) 

A different choice of the standardization of the wavefunction for \z\ — )• oo, the 
standing wave boundary conditions, will be more convenient for numerical calculations. 
In order to make explicit use of the invariance of the problem under the coordinate 
inversion r — )■ —f we first define new sets of reference functions in the z motion with 
well defined behavior under the z — ?■ ~z refiection vr^. Namely, diagonal matrices of 
functions F'^^ and G'^^ regular and irregular at the origin are defined for open channels 
and parity vr^ = ± by 

= kn^^"^ cos{knz) I G:^{z) = kn^^"^ sm{kn\z\) 

(z) = kn^^"^ sm{knz) ' I G;;{z) = -kn^^'^sign{z)cos{knz). 
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The corresponding reference functions for closed channels are 

F+{z) = cosh{knz) I Gf^{z) = - exp{-K\z\) ^^^^ 

P- (z) = smh.{knz) ' I G;'^{z) = -sign{z)exp{-kn\z\). 

With these definitions, total parity of product functions (prikF^ and (pnAGt 
p = 7r2(— 1)'^. We can now define a real reactance matrix K'^"^ for each vr^ = ± and 
angular momentum A according to the asymptotic behaviour [21j 

0nA (P, ^) F^iz) - J2 K'n^n'K (p, ^) Gt{z) , \z\ ^ CX).(ll) 

n' 

If the transition matrix is also decomposed into symmetric and antisymmetric 
components P"', 

(sign(z)fc„, ^ = + sign(^)/„-„(E) (12) 

the matrix relations 

=iK'^- (I-iK^-)-^ (13) 

obtained comparing ([71) and ( ITTi) can then be used to determine the in terms of K'^^. 
Finally, we will make use below of the unitary and symmetric scattering matrix, defined 
for the two parities as S''^ = (I + i^'^') (I - zK''-)"\ 

2.2. Wigner laws 

The scattering quantities obey at low energy Wigner laws, valid when the asymptotic 
de Broglie wavelength is much larger than the range of the effective potential in a given 
channel. Following Ref. [22], a collision channel in which the collision energy tends to 
zero will be termed new. Similarly to the 3D case, in which Wigner laws depend on 
the angular momentum quantum number, even in the present ID configuration it is 
convenient to treat separately the tt^ = ± parity cases. An analysis along the lines of 
Ref. [22j gives the asymptotic behaviour, 

/+„~-l + 0(fc„) , /-„~A;„. (14) 

implying that for vanishing collision energy at the ?7,-th collision threshold a Landau state 
will be reflected with unit probability, irrespective of the value of A. Moreover, the odd 
transition matrix element becomes negligible with respect to the even component , such 
that particles will be scattered with equal probability in forward and backward 

directions. 

Inelastic processes consist in a transition to a Landau state of different transverse 
energy, and the corresponding scattering amplitudes obey for both parities the 
asymptotic law 

fn'n ~ (15) 

if n is new. The transition probability to a newly energetically accessible channel is 
then vanishing at threshold both in reflection and transmission. A similar relation 

fn'n - VK' (16) 
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holds in the case where n' is new. 

The Wigner laws allows one to define for positive parity an effective coupling 
constant giT>{kn) finite at threshold — )■ 0), which embeds the details of the short 
range dynamics in a single parameter 

fnn ~ ^27 (-'-'^) 
, n Kn 
1 — I 

Such gmikn) can be used for instance to parametrize a ID pseudopotential Vp, i.e. 
a zero-range interaction Vp = gi]:,S{z) which generates the same elastic scattering 
amplitude as the full problem at momentum kn. An associated one dimensional 
scattering length can be defined for one open channel as the limit value oid = 
— limic^Q f] / k , where ri{k) is the scattering phase shift extracted from the asymptotic 
wave function \E'^"*"(r) oc sin(/c|z| + ?7)0oo (p) [23]. Its relation with the scattering 
amplitude is 

aiD huai^o /.qn 
= . [IH) 

O-ho giD 

For a purely elastic collision guy is real, otherwise its non-vanishing imaginary part 
accounts for inelastic loss processes. Note that unlike the three-dimensional case, 
where is proportional to g^^, here the scattering length is inversely proportional 
to the coupling constant. We will restrict most of our discussions to giD, since this is 
the quantity usually measured in experiments [21]. Moreover, pseudopotentials with 
coupling strength gn) are often used to model two-body interactions in many-body 
theories. A numerically exact approach for the determination of the scattering amplitude 
and hence of the coupling constant is the subject of section [3l 

2.3. Feshbach resonances in quasi ID 

Feshbach resonance theory is well established for three dimensional collisions. We shortly 
recall here its main elements in order to point out some differences with respect to 
the present quasi one- dimensional configuration. Let us focus for definiteness on the 
case of one open channel and one isolated resonance. We will consider positive parity 
scattering and drop the tt^ parity label for notational convenience. A formal separation 
of background and resonant contributions to the scattering matrix can be performed 
following the standard Fano- Feshbach approach [25| [26] 

2{Ej-E^-Ce)^zje .... 

where S'bg is the off-resonance scattering matrix, Eq the energy location of the molecular 
state responsible for the resonance, ( and 7 the resonance shift and width, respectively. 
Near resonance the quantity Eq typically presents a smooth linear dependence on the 
external parameters of the problem, the magnetic field or the potential depth in the 
present case. 
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At low momentum, according to the Wigner laws S ~ (—1 — 2ikaiD), S'bg ~ 
(— 1 — 2z/caiD,bg)5 Ce ~ const and je ~ k. Taking the -E — ?■ limit, we obtain the 
resonantly modified oid as 

aiD = aiD.bg (^1 - '^^^r^ ' 

where we have defined A = limfc_j.o 7_E/(2A;aiD,bg)- Note that on resonance {Eq = —(q) 
the ID scattering length diverges and according to f llSp the system becomes non 
interacting, qid — )■ and the transmission coefficient tends to one. This behaviour 
should be contrasted with the infinitely strong interaction obtained in three dimensions 
as the molecular state crosses the collision threshold. Note however that near resonance 
aij) will cross zero, and Qid will diverge, if Eq + (q = A. With this remark, in the 
following we will refer to the divergence of the effective coupling constant as to a zero- 
energy resonance situation. 

3. Numerical methods 

We solve our problem in the framework of the time independent close coupling approach 
to potential scattering. In this method the full wavefunction is expanded in a suitable 
complete basis depending on a reduced set of coordinates. The choice of the basis 
function is inspired by the physics of the problem. Let us define a spherical system of 
coordinates {r, 6, (f} with origin on the collision center and polar axis in the z direction. 
We identify an "inner region" r < Rq where the isotropic interatomic interaction initially 
dominates. The criterion for choosing the boundary Rq is given below. In this internal 
region an expansion in spherical harmonics is used 

^fl(^1= E'^"'/^"W^^A(r), (21) 

^>|A| 

where the sum runs only over even or odd values of £ for positive or negative parity, 
respectively. The coupled channel Schrodinger equation takes the form 



-I— -7^ + V + C - Awl 



f(r) = M(r). (22) 



v,{r) = ZZ2 + ^(^) (23) 



The effective centrifugal potential v has a simple diagonal form with elements 

h ^i{i + 1) 

2mr^ 

whereas the coupling matrix C is formed by the matrix elements of the transverse 
oscillator potential 

2 2 

Tf) 111 T 

= -^{Y,a\ sm\e)\YeA). (24) 

Eq. fl22l) is solved up to Rq using the Johnson-Manolopoulos propagator [27], which 
computes the log derivative matrix Z = (^f)f~^ along a reaction coordinate using 
embedding propagators. 
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In the external r > Rq region it is more natural to use cylindrical coordinates and 
an orthogonal basis of Landau oscillator states, Eq. (jl]) 

^nii^ =Y.K>n{z)<Pn'K{p.V)- (25) 

n' 

The coupled channel Schrodinger equation in the cylindrical basis becomes 



2m dz"^ 



Er +U 



h(^) = Eh{z). (26) 



The diagonal matrix E^, is composed by the Landau levels ([3]) and coupling u arises 
from the spherical atom-ion interaction V . In practice we will always switch to the 
cylindrical representation at distances where the V{r) = —Cir'"^ asymptotic form is 
accurate. Moreover, in the same region the size of transverse oscillator states is much 
smaller than the radial coordinate r, such that a series expansion in the small parameter 
(^ho/^ can be performed. We find matrix elements 

Un'n = {(pn'A\V\(f)nA) = ^1 - ^ (0n'A | ^wa) H ^ • (27) 

The second term in the expansion is rapidly decreasing ~ z~^ and in the present 
calculations we only retain lower order ~ z~'^ diagonal coupling elements Unn for both 
open and closed channels. 

The representations (!?!]) and are matched on the spherical surface 5*0 defined by 
r = Rq. We find that to insure optimal numerical convergence one should choose Rq in 
such a way that the relevant Landau states are localized in the two hemispheres without 
significant overlap on the equatorial plane 0. The matching procedure is analogous to the 
one used in reactive scattering [28l [29] . We anticipate that the change from spherical to 
cylindrical geometry is at the origin of the peculiar scattering and resonance phenomena 
described in the next section. 

In the external region r > Rq we integrate inward the equations ( l26l) from very 
large values of z, with starting value the asymptotic boundary condition (!9|)- (fT0|) . Let 
us denote the corresponding matrix solutions and G^, which retain diagonal form 
due to the neglect of the off diagonal elements in the coupling u. The scattering solution 
is represented as 

^nli^ = T.[^^('^)^-'n-G^,iz)KX]<l>n'A{p,^) , r>Ro. (28) 

n' 

Continuity of the scattering wavefunction and of its normal derivative on the 
spherical surface 5*0 is used to determine K in terms of the logarithmic derivative matrix 
Z. The projection matrices 

J^in = {yiA\Fn4>nA) So 

J-;,= (F,a||:(F„0„a))5o (29) 

I Typically, Rq — lOoho insures a good convergence for collisions starting in the lowest Landau levels. 
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needed to this aim |29j are numerically evaluated using Gauss-Legendre quadratures. 
Similar projection matrices Q^n and involving the irregular reference functions are 
defined by replacing F with G on the rhs of (l29i) . The reactance matrix is finally 
obtained for each parity using the equation [29] 

¥. = {zg -g'y\z7 -7'). (30) 

4. Results and discussion 

In this section we apply our formalism to atom-ion diffusion with A = 0, a somewhat 
idealized situation to test our theoretical formalism which can however be in principle 
experimentally realized. For instance, instantaneous switching of a magnetic field in the 
propagation direction of an ion beam would project the plane wave describing the beam 
on a superposition of Landau states with different n and A = without affecting the 
longitudinal momentum k. A successive measurement of the transverse energy could be 
used to select the desired initial asymptotic state n for the collision. 

The default molecular parameters we choose for our model correspond to the 
Na+Ca"*" system. However, previous work has shown that the short range molecular 
physics can be accurately described in the ultracold regime by a single quantum defect 
parameter independent of energy and angular momentum, which can be identified with 
the s-wave three-dimensional scattering length pQ]. Within this approximation, our 
results will be cast in an essentially system independent form by presenting them in 
terms of the only independent dimensionless groups 

Cts E Oho 

that can be built using the physical parameters of the problem. Note that in this work 
we set m equal to the reduced mass of the system, for future comparison with the 
free-space results |10j . 

We first consider collisions starting in the lowest Landau level n = aX zero 
collision energy and calculate the effective coupling constant gmiko — J- 0). We make 
the three dimensional scattering length vary from infinite negative to positive infinite 
values by artificially modifying the interaction potential. The variation of in units of 
ajio will be mapped to a finite interval by introducing the normalized quantum defect 
^ = (2/7r) arctan(as/aho)- This analysis will be repeated for three values of a^o. We will 
then analyze the behaviour with the collision energy of the scattering observables for 
some selected values of a^o and a^. 

4.1. Weak field regime 

Let us first fix a value of B such that the size of the lowest Landau state is significantly 
larger than the interaction range. We will refer to this condition as the weak field 
limit, and choose more specifically a^o = lOR* {B = 10.9 G). A salient feature observed 
in Fig. [1] is the strong divergence of the effective coupling constant at a characteristic 
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value of as- That is, the effective ID couphng strength can be tuned to virtually any 
value by modifying the value of Og with respect to aho, for instance through a variation 
of the applied magnetic field. We stress that this zero energy resonance arises from 
the quantization of the continuum in Landau levels and does not depend on resonance 
phenomena preexisting in free space. As such it can be termed a geometric resonance 
effect. 



1 — I — I — I — I — I — I — I — I — I — I — r 




J I I I I I I I L 



Figure 1. Numerically computed effective coupling constant for Ca++Na collisions 
in a weak magnetic field with = fOi?* as a function of the normalized scattering 
phase ^ = (2/7r) arctan(as/aho) for axial angular momentum A'^^ = 0+ (full line). 
Vertical dotted lines denote the position where giD diverges (aiD = 0). Energy 
independent end energy dependent pseudopotential approximations (see text) labeled 
I and II respectively are also shown. 



Interestingly, since in the weak field limit the range of the potential can be neglected 
with respect to the transverse extension of the Landau state, the interaction V can 
be expected to be well described by a zero range 3D pseudopotential. With this 
approximation, Eq. (!26|l for A = reduces to an analytically solvable model [231 EH] • 
The effective coupling constant of Ref. [23] is 



2- 



, - .l + ^C(l/2) , (32) 

where C, is the Riemann zeta function and C(l/2) = —1.46035. The zero-order result 
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(7iD = 2ashuj, obtained by neglecting the second term in the parenthesis on the rhs, 
simply represents the renormalization of Og obtained in the Born approximation by 
averaging the three-dimensional pseudopotential over the initial transverse Landau 
state. However, Eq. ( !32|) also shows that a confinement- induced resonance arising 
from virtual excitations to asymptotically closed Landau states should appear for 
as = -aho/C(l/2) ^ 0.6848 a^,. 

This nontrivial prediction is confirmed by the comparison with our numerical 
calculation shown in Fig. [H The small shift in the resonance position can be understood 
by realizing that the result (132|) is exact only in the limit of vanishingly small axial 
collision energy. However, when the atom and the ion are close together the relative 
energy E in their unconfined 3D motion does not vanish but is rather equal to the zero 
point energy hw, see Eq. ([3]). This correction can be introduced by using an energy 
dependent 3D peudopotential [31], which replaces the zero energy parameter with its 
finite energy counterpart — tan((53D)//c3D- Here S^b is the s-wave phase shift at energy 
E and k^jy = a/ 2mE / is the wave vector in the three-dimensional motion. With 
this modification our numerical calculation becomes virtually identical to Eq. f l32|) thus 
confirming the accuracy of the pseudopotential approximation even for a long-range r~'^ 
potential under the appropriate conditions. 




Figure 2. Same as figure [1] but in the intermediate field case Oho = R* ■ The result of 
a purely s-wave calculation with the full molecular potential is also shown (see text) . 
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Since the zero point energy is small with respect to the centrifugal barriers in 
Eq. fl22|) one can expect, as indeed proven by the agreement with the s wave model 
( 132|) . that diffusion in partial waves £ > should be significantly suppressed. Our 
numerical approach allows us to check this property by a direct calculation. To this 
aim, we set = in the f£>o components of the interaction matrix v. The non s-wave 
potentials are also turned off in the external region by using in the projection matrix 
elements f l29|) with £ > the unperturbed waves F and G in the place of F and G. Note 
that in spite of the fact that dynamically a single partial wave contributes to scattering, 
one cannot simply use a unique £ = in the expansion ( 12T|) , since a sufficient number of 
partial waves is in any instance cinematically necessary in order to represent the relevant 
Landau states on the spherical surface 5*0 of radius Rq. The result of this purely s-wave 
calculation is not included in Fig. [T] as it is virtually indistinguishable from the exact 
result. 

4-2. Intermediate field regime 

Let us now consider the variation of qid with in the intermediate field regime aho = R* 
[B = 1.09 X 10^ G). Like in weak field conditions, a single geometric resonance is observed 
in the numerical simulation, see Fig. [2l The pseudopotential result in Eq. ( l32l) in 
the energy-independent (I) and energy- dependent (II) formulations is only qualitatively 
correct. The purely s-wave calculation shows a significant shift in the resonance position 
from both the exact calculation and the pseudopotential models. 

These findings can be rationalized as follows. The discrepancy between the 
pseudopotential and the purely s-wave result demonstrates the breakdown of the 
assumption of a zero-range potential, that is the importance of the long-range interplay 
between atom-ion and magnetic field potentials. The disagreement between the s-wave 
and the exact result is due to the important contribution of higher order partial waves. 
This contribution does not vanish even in the — )■ limit because of the zero point 
energy. Its value (~ 3/iK), a factor 10^ larger than in the weak field case, and the long 
range nature of the potential insure that a few partial waves are scattered from the 
potential. More quantitatively, by artificially varying the number of partial waves as 
described above we find that a calculation restricted to £ = — 4 determines with high 
accuracy gn). The good agreement in the resonance position between the exact result 
and the pseudopotential (I) approximation should be regarded as merely coincidential, 
stemming from the compensation of long-range and higher order partial waves effects. 

4-3. Strong field regime 

Finally, Fig. |3] shows the behaviour of gin for a strong field case Oho = O.IR*, a tight 
confinement condition in which the potential range is larger than the Landau orbit 
radius. The required magnetic field intensity is of 5 = 1.09 x 10^ G, a large yet 
experimentally attainable value. In this conditions both the pseudopotential and the 
purely s-wave model completely fail to reproduce the numerical results. Several zero- 
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Figure 3. Same as figure [T] but in the strong field case aho = O.li?*. The 
pseudopotential approximation breaks down in this situation and is not shown. Note 
that for graphical purposes the normalized quantum defect on the horizontal axis is 
defined in terms of R* rather than of Oho- 



energy resonances are now available as (or B) is made to vary. The nature of such 
resonances can be understood from test calculations in which the atom-ion potential is 
suppressed in selected partial waves. Such calculations show that each broad feature in 
the QiD of Fig. |3] arises from the combined effect of a relatively large number of partial 
waves £ = — 12. That is, due to the relatively large zero point energy {hw ^ 300/iK) 
the wave function picks up a significant phase shift at short range from different partial 
wave effective potentials. The narrowest features observed near ^ = — 1 can in contrast 
be assigned to a single or few i > 4 angular momentum values, suggesting that most of 
the amplitude of the resonant bound state is localized in the potential well region. An 
analysis of the bound states of the systems will be needed to assign more precisely the 
features observed in the present scattering calculations. 

4.4- Behaviour with collision energy 

The zero energy resonances described so far may shift as a function of the potential 
strength into the scattering continuum giving rise to possible resonance effects 




Figure 4. Elastic reflection probability for an ion in the lowest Landau level as a 
function of collision energy in the intermediate field case aho — R*- Axial angular 
momentum is A = and typical values of the tridimensional scattering length 
as = ±i?* are selected. Odd and even contributions to the reflection probability 
are separately shown for the as = R* case. Resonance and threshold phenomena in 
the reflection probability are discussed in the text. 



observable as a function of collision energy. Let us consider in detail two potentials 
supporting "typical" scattering lengths = ±_R* in the intermediate field case Oho = -R*. 
Fig. m shows the elastic reflection probability Rqq for a scattering event starting with the 
ion in the lowest Landau state 0oO) as a function of collision energy. A decomposition 
of the scattering amplitude in even and odd components shows that the small value of 
the reflection coefficient i?oo = l/oo ~ /ooP between the n = and 1 thresholds is due 
to interference of parity amplitudes of relatively large magnitude, see Fig. HI A similar 
situation has been reported for neutral atom scattering from short-range potentials in 
[32]. Fig. [5] stresses the fact that the Wigner regime near the n = threshold at which 
RoQ — )■ 1 is correctly reproduced by the numerical calculation but is extremely narrow, 
for as = —R* in particular. 

For as = R* a. first peak at which i?oo reaches one is observed below the opening of 
the n = 1 level with energy fiWc- This feature arises from the positive parity part f+ of 
the scattering amplitude, and can be supposed to be a resonance arising from trapping 
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Figure 5. Near threshold elastic reflection probability for as 
parameters are as in Fig. HI 
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±R* . Systems 



To confirm the resonant character of this feature we calculate the time delay matrix 
introduced by Smith [33], defined for each parity as 



dE 



(33) 



Note that Q is expressed in terms of the open channel scattering matrix, such that 
its dimension equals the number of energetically accessible channels. The evolution 
with energy of the largest eigenvalue q of the lifetime matrix can provide in particular 
information on the lifetime of the complex formed during a collision. For instance, a 
narrow isolated resonance far from a collision threshold is characterized by a time delay 
following a familiar Lorentzian profile [331 El] • In general, at very low collision energy 
threshold effects produce deviations from this simple picture making the analysis less 
direct [35]. In the present case the situation is simpler, since the lifetime obtained from 
Q"*" shows a pronounced bell shape profile of width much smaller than its central energy. 

This analysis provides conclusive evidence of the resonant nature of the feature 
observed in the reflection probability (see Fig. [6]) and allows us to assign positive parity 
to the associated metastable state. As already noted in [35] the time delay diverges 
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positively or negatively at zero energy with a 2maiY)/{hk) law, corresponding to the 
classical time to span the 2ai£) distance. 



300 
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Figure 6. The two largest eigenvalues of the time delay matrix as a function of 
collision energy for the = R* collision in the intermediate field case of Fig. 2] 
Only 71^ = + eigenvalues are shown. The long-lived isolated peaks are conclusively of 
resonant origin. 



Below the second n = \ and third n = 2 thresholds additional rapid variations with 
energy of smaller amplitude appear in i?oo- Again, the time delay analysis shown in 
Fig. |6] conclusively confirms the resonant nature of these features. The dispersive-like 
shape observed in the reflection probability arises from interference with background 
scattering, whereas the small modulation amplitude is explained by quenching due to 
n = — )■ ra' > inelastic processes. 

Quite strikingly, such resonances seem to occur at similar locations with respect to 
each threshold. This property can be qualitatively interpreted using the pseudopotential 
analytical method of Ref. [36]. In this approach the resonance below the first 
threshold is assigned to a molecular state generated by the zero-range interaction in 
the n > subspace of closed channels. By a straightforward generalization consisting 
in replacing the shift operator = J2n\^n+i,o) {(pn,o\ below Eq. (10) of |36j by 
= ^„ |0n+m,o)(0n,o|, oue cau show that if a resonant state is present below the 
n = 1 level there has to be a corresponding one below the mth threshold. Similar 
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Figure 7. Inelastic reflection and transmission probability for an ion in the lowest 
Landau level as a function of collision energy in the intermediate field case aho = R* 
for A = and = -R*. 

regularities in the Landau resonance spectrum have also been proved in low-energy 
electron photodetachment using a simplified model which assumes resonant states 
confined in the z = plane [37] . 

The Qs = —R* refiection coefficient shows an evident cusp at the opening of the first 
threshold. Near the higher thresholds, rounded step features are also evident. These 
features arise from the non analyticity of the scattering amplitude near the opening of 
a new threshold and do not have resonant origin. In fact, as shown by Newton both 
sharp and rounded features can appear at threshold openings |38j . 

We also show in Fig. [7] inelastic transition probabilities from the lowest to excited 
Landau states for = —R*. Note that T^o and Rno vanish at the nth threshold, in 
agreement with the Wigner laws, and are otherwise of significant magnitude. The T„o 
and Rno differ in general, showing once again that scattering from odd partial waves 
is non negligible in the intermediate ffeld case, and introduces a signiffcant forward- 
backward asymmetry. As expected the i?„o and T„o do not vanish at thresholds m 7^ n, 
but present cusps and rounded step singularities similar to the one observed in the 
elastic probabilities. 

Finally, we compute the elastic reffection probability in the strong ffeld limit 
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Figure 8. Elastic reflection probability for an ion in the lowest Landau state as 
a function of collision energy in the strong field case Oho — O.li?*. Axial angular 
momentum is A = and typical values of the tridimensional scattering length 
as = ±R* are selected. 

(Fig. |8]). The rich resonance structure characteristic of zero energy scattering has a 
finite energy counterpart. In fact, multiple resonance peaks are now observed between 
Landau thresholds. The narrowest features are most likely associated to bound states 
of relatively large (. components. At variance with the intermediate-field case, a 
decomposition of the transition matrix identifies the resonance peaks as originating 
from both positive and negative parities. Cusp and rounded steps stemming from the 
non analytical behaviour of the scattering matrix and from interference effects between 
parity components of the scattering amplitude are also visible near the opening of new 
energy thresholds. 

5. Conclusions and outlook 

In conclusion, we have presented a rigorous computational approach to study atom- 
ion collisions in the presence of a magnetic field. The weak-field case maps to a well 
understood analytical model. In contrast, we have shown that a novel regime of strong 
confinement can be reached in atom-ion systems at experimentally attainable values 
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of the magnetic field. A series of zero- and finite-energy resonances associated to higli 
order partial waves appears in this case. The resonance phenomena we have analyzed 
in this work should play an important role in applications based on magnetic control of 
the atom-ion dynamics. 

In a forthcoming paper we will present our results for the qualitatively different 
case of |A| > collisions and consider the case of a scattering center not necessarily 
placed on the axis of the cyclotron orbit. In perspective, it will be interesting to 
investigate the interplay of pre-existing atom-ion magnetic Feshbach resonances and of 
Landau quantization. Moreover, in cases where the atom cannot be treated as infinitely 
massive, quasi separability of the center of mass motion might give rise to still unexplored 
resonance phenomena. 
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